Right-Handed Neutrinos as the Dark Radiation: Status and Forecasts for the LHC 
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Precision data from cosmology (probing the CMB decoupling epoch) and light-element abundances 
(probing the BBN epoch) have hinted at the presence of extra relativistic degrees of freedom, the 
so-called "dark radiation." We present a model independent study to account for the dark radiation 
by means of the right-handed partners of the three, left-handed, standard model neutrinos. We 
show that milli-weak interactions of these Dirac states (through their coupling to a TeV-scale Z' 
gauge boson) may allow the vr's to decouple much earlier, at a higher temperature, than their 
left-handed counterparts. If the vr's decouple during the quark-hadron crossover transition, they 
are considerably cooler than the vl 's and contribute less than 3 extra "equivalent neutrinos" to the 
early Universe energy density. For decoupling in this transition region, the 3 vr generate AN V = 
3(Tv R /T VL ) 4 < 3, extra relativistic degrees of freedom at BBN and at the CMB epochs. Consistency 
with present constraints on dark radiation permits us to identify the allowed region in the parameter 
space of Z' masses and couplings. Remarkably, the allowed region is within the range of discovery 
of LHC14. 



Big-bang nucleosynthesis (BBN) is remarkably succes- 
full in predicting the relative abundance of light elements 
as a function of two fundamental parameters: the baryon 
density of the universe, Ob/i 2 , and the number of "equiv- 
alent" light neutrino species, iV c ff [HQ In fact, until re- 
cently BBN provided the only constraint on these pa- 
rameters. Discovery of primordial anisotropics in the 
cosmic microwave background (CMB) has granted a su- 
perior test-bed for precision constraints on fundamental 
parameters in cosmology. This powerful test-bed can be 
used to assess whether new physics model predictions are 
simultaneously consistent with BBN and CMB observa- 
tions. Of interest here is the capacity toprobe right- 
handed neutrino milli-weak interactions 13j-|7|, which are 
predicted in various extensions of the standard model of 
particle physics. 

In this Letter we construct a model independent tem- 
plate for placing upper and lower bounds on the mass 
of an extra Z' gauge boson, which allows for milli-weak 
interactions of the right-handed partner of the Dirac neu- 
trino. A critical input for such an analysis is the relation 
between the relativistic degrees of freedom (r.d.o.f.) and 
the temperature of the primordial plasma. This rela- 
tion is complicated because the temperature which is of 
interest for right-handed neutrino decoupling from the 
heat bath may lay in the vicinity of the quark-hadron 



1 The discussion here is in the context of the usual concordance 
cosmology of a flat universe dominated by a cosmological con- 
stant, with Qyv ~ 0.7 and a cold dark matter plus baryon compo- 
nent Q m ~ 0.3; the Hubble parameter as a function of redshift is 
given by H 2 (z) = Hi [f2 m (l + z) 3 + QaL normalized to its value 
today, H ~ 100 h km s" 1 Mpc -1 , with h ~ 0.71 0|. 



cross-over transition. In a previous publication Q, use 
was made of a detailed lattice study to connect the 
temperature to an effective number of degrees of free- 
dom. Very recently, one of us has provided an analysis 
in which the decoupling of the extra relativistic degrees 
of freedom may occur well beyond the cross-over tem- 
perature [1, Q. In that case a general connection be- 
tween the effective number of degrees of freedom and the 
right-handed neutrino decoupling temperature is needed. 
To this end, we employ the results of Ref. to find 
the post-e 1 * 1 annihilation ratio of the temperatures of the 
right-handed and left-handed neutrinos, T VR /T VL , which 
is then used to predict the enhancement to the effec- 
tive number of degrees of freedom in the early Universe, 
AN V = Z{T V JT VL Y < 3. 

The formulation presented here allows for an immedi- 
ate test for the potential of any model to account for any 
extra neutrino degrees of freedom. For illustration, we 
analyze several candidate models. Before proceeding we 
provide a brief and concise overview of the current ob- 
servational constraints on the number of light neutrino 
species. 

Over the past few years, the Wilkinson Microwave 
Anisotropy Probe (WMAP) the Atacama Cosmol- 
ogy Telescope (ACT) [H, and the South Pole Telescope 
(SPT) [l3| have each provided evidence for a "dark" rel- 
ativistic background (a.k.a. dark radiation). Parame- 
terized in terms of the number of relativistic degrees of 
freedom the data seem to favor the existence of roughly 
one extra effective neutrino species. Specifically, the pa- 
rameter constraint from the combination of WMAP 7- 
year data, the latest distance measurements from the 
baryon acoustic oscillations (BAO) in the distribution of 
galaxies 14 1, and precise measurements of H [l{| lead 
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FIG. 1. Comparing the BBN predictions of N e s and Qsh 2 
with those from various CMB determinations: BBN D + 4 He 
(red filled triangle) [1, BBN D + WMAP7 [U (red open 
triangle), WMAP7 [n[_[blue filled square), ACT (green 
filledpentagon), SPT [13] (purple filled circle), SPT + Clus- 
ters [Ts| (purple open circle). Taken from Ref. 



to N cS 



4.341° ;fg (68%CL) 



II II ] . Similarly, a combi- 



nation of BAO and Hq with data from the ACT yields 
N cS = 4.56 ± 0.75 (68%CL) 0, whereas data collected 
with the SPT combined with BAO and Ho arrive at 
N cS = 3.86 ± 0.42 (la) 0. 

Turning now to the BBN determinations, we note that 
the observationally-inferrcd primordial fractions of bary- 
onic mass in 4 Hc had been favoring N e g < 3 16 1. Un- 



expectedly, two recent independent studies determined 
a larger 4 He fraction and the updated effective number 
of light neutrino species is reported as N e g = 3.80+Q'yQ 

(2a) nam. 

Several recent papers have presented results from com- 
bined analysis of the BBN and /or CMB data, but based 
on different priors (see e.g. |19l423j ). Throughout we 
adopt N e s = 3.7ltQ 45 (la), which represents a conser- 
vative choice based on an extensive BBN analysis @ . As 
shown in Ref. [1] and display here in Fig. [1] this analysis 
is in agreement with various CMB observations. In the 
standard electroweak theory, N B g — 3.046 (the difference 
from 3 being mainly due to partial heating of vl by e + e~ 
annihilation (24[), yielding 



AN„ = 0M± o H (la) 



(1) 



Several explanations have been proposed to explain 
such a possible AN V excess. These include: (i) mod- 
els in which the extra relativistic degrees of freedom are 



related to possible dark matter candidates [25j-|30(; (ii) 
models based on active-sterile mixing of neutrinos in a 
heat bath pjlL l32|: (Hi) models based on milli-wcak in- 
teractions of right-handed partners of three Dirac neu- 
trinos Q ■ In this work we confine our discussion to case 
(Hi). 

We begin by first establishing, in a model independent 
manner, the range of decoupling temperatures implied by 
the BBN or CMB observations. The effective number of 
neutrino species contributing to r.d.o.f. can be written as 
N e g = 3[1 + (T VR /T VL ) 4 ] ; therefore, taking into account 
the isentropic heating of the rest of the plasma between 



T^ c and T^ c decoupling temperatures we obtain 



4/3 



(2) 



where g(T) is the effective number of interacting (ther- 
mally coupled) r.d.o.f. at temperature T; for example, 
g(T^ c ) = 43/4 (Hi For the particle content of the 
standard model, there is a maximum of g{T^ c ) = 427/4 
(with T^° c > m t0 p), which corresponds to a minimum 
value of AN V = 0.14. For the subsequent study, we 
adopt the determination of g(T) given in Q based on 
the results of [l(|. Then using Eq. © we obtain a re- 
lation between AN V vs. X^ c , which is shown in Fig. [2] 
From this curve we determine the range of decoupling 
temperature: T^ c = 0.174±J;g|g GeV. 

The physics of interest then will be taking place at en- 
ergies in the region of the quark-hadron crossover tran- 
sition, so that we will restrict ourselves to the follow- 
ing fermionic fields, and their contribution to r.d.o.f.: 
[3u R ] + [3d R ] + [3s R ] + [?>i>l + e L + fi L ] + [e R + fi R ] + 
[3ul + 3o?l + 3s l] + [3i/fl]. This amounts to 28 Weyl 
fields, translating to 56 fermionic r.d.o. fH 

The right-handed neutrino decouples from the plasma 
when its mean free path becomes greater than the Hubble 
radius at that time. To determine T^ c , wc first calculate 
the vr interaction rate 



r(T) 



( 9 
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i— 

where A/i is the number of chiral states, 

1/4 
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(3) 



(4) 



If relativistic particles are present that have decoupled from the 
photons, it is necessary to distinguish between two kinds of g: g p 
which is associated with the total energy density, and g B which 
is associated with the total entropy density. For our calculations 
we use g = g p = g a . 

In principle, the contributions to the vr interaction rate T(T) 
from the c quark and r lepton should be included for the lowest 
value of A7V„, which in this paper corresponds to a decoupling 
temperature of 1.5 GeV. However, one can easily verify that this 
inclusion will not be visible in Fig. [3] 
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FIG. 2. Relation between AiV^ vs. T° 



gi are the chiral couplings of the Z' gauge boson, and the 
cosntant K. = 0.5 (2.5) for annihilation (annihilation + 
scattering) 0. Armed with (|3]) we determine T^° c via 
the prescription 



r(T dcc } 



where 



Mpi I Ai\T„ 



3/8 



(5) 



(6) 



Substituting © and ([5]) into (|S|) we obtain 



3/32 / 13.28 



1/4 



Mpi AC (T„ d -) 3 



(7) 



For a given value of AA„, Q conveniently yields a 
straight line plot ofgvs. Mz> ■ In Fig.[3]we provide graphs 
corresponding to the central value and la limits for the 
values of A A„ given in (TTJ) . The hatched region between 
the highest and lowest lines reprersents the ~g — Mz> pa- 
rameter space consistent with the AN V measurement. 

To illustrate we calculate ~g for two candidate models. 
The first is a set of variations on D-brane constructions 
which do not have coupling constant unification. The 
second are two U(l) models (U(l)^ and U(l)x) which 
are embedded in a grand unified exceptional Eq group, 
with breaking pattern 

E 6 -> 50(10) x U(l)^ -> 5*7(5) x U(l)^ x U(l) x . (8) 



The latter two are interesting because they provided a 
test basis for Z' searches at ATLAS and CMS (35|. 
For each of the Eq models we may write gi in Q as 
9i = .9o0ij_where in conformity with grand unification 
we follow \& and choose 



9o 



gi t&nOw ~ 0.46 , 



(9) 



with (72 the SU{2)l coupling. The charges Qi for the 
different fermions are conveniently tabulated in 0] . 

In the D-brane construction, the Weyl fermions live at 
the brane intersections of a particular 4-stack quiver con- 
figuration: U(3) c x SU(2) L x U(l) Ift x U{1) L The 
resulting U(l) content gauges the baryon number B [with 
U(1)b C U(3)c], the lepton number L, and a third addi- 
tional abelian charge Ir which acts as the third isospin 
component of an SU (2)r. Contact with gauge structures 
at TeV energies is achieved by a field rotation to couple 
diagonally to hypercharge Y. Two of the Eulcr angles are 
determined by this rotation and the third one is chosen 
so that one of the U(l) gauge bosons couples only to an 
anomaly free linear combination of In and B — L [37| . 
Of the three original abelian couplings, the baryon num- 
ber coupling is fixed to be y/1/6 of the QCD coupling 
at the string scale. The orthogonal nature of the rota- 
tion imposes one additional constraint on the remaining 
couplings [38( . Since one of the two extra gauge bosons 
is coupled to an anomalous current, its mass is 0(M S ), 
as generated through some Stiickclbcrg mechanism. The 
other gauge boson is coupled to an anomaly free current 
and therefore (under certain topological conditions) it 
can remain massless and grow a TeV-scalc mass through 
ordinary Higgs mechanisms [39l |. We consider two ex- 
treme possibilities in which the TeV-scale Z' gauge bo- 
son is mostly or mostly B — L. The chiral couplings 
{gi) of these gauge boson are tabulated in 37j. We also 
consider a D-brane construct with TeV-scale string com- 
pactification (the chiral couplings of this model are given 



in Table IV of [38j]). Details of these assignments are 
given in the figure caption. Termination of the lines on 
the left reflects the LHC experimental limits on the mass 
of the gauge boson, using null si gna ls for enhancements 
in dilepton 13, 35| and dijet 0, 41 1 searches. 

The cosmology results from the Planck satellite would 
allow determination of N e g with a standard deviation 
of about 0.2 42, Hijl, whereas the future Large Synoptic 
Survey Telescope (LSST) could determine N c g with a 
standard deviation of about 0.1 [44[. With this enhanced 
sensitivity the hatched region will collapse to a line and 
intersect for any given model its horizontal curve at the 
mass of the Z'. 

Note added: After this work was finished a paper ap- 
peared on the arXiv with a comprehensive study on dark 
radiation of Eq models [45[. Our results are completely 
consistent with those of Ref . 45 1 . 
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FIG. 3. The green cross-hatched areas show the region allowed from decoupling requirements to accommodate BBN and CMB 
eras. Each of the horizontal lines refers to a particular model: from the top the E& Z' x , a D-brane model in which Z' is mostly 
B — L, a D-brane model in which Z' couples mostly to the third component of a right-handed isospin, a D-brane model with 
TeV-scale strings, the Ee ZL. Termination of the lines on the left reflects the LHC experimental limits on the mass of the gauge 
boson. The left and right figures show the condition on decoupling for loss of chemical and thermal equilibrium, respectively. 
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